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ABSTRACT 

No optical afterglow was found for the dark burst GRB 981226 and hence no absorption redshift 
has been obtained. We here use ground-based and space imaging observations to analyse the spec- 
tral energy distribution (SED) of the host galaxy. By comparison with synthetic template spectra 
we determine the photometric redshift of the GRB 981226 host to be Zphot = 1.11±0.06 (68% confi- 
dence level). While the age-metallicity degeneracy for the host SED complicates the determination 
of accurate ages, metallicity, and extinction, the photometric redshift is robust. The inferred Zphot 
value is also robust compared to a Bayesian redshift estimator which gives Zphot = 0.94 ± 0.13. The 
characteristics for this host are similar to other GRB hosts previously examined. Available low reso- 
lution spectra show no emission lines at the expected wavelengths. The photometric redshift estimate 
indicates an isotropic energy release consistent with the Amati relation for this GRB which had a 
spectrum characteristic of an X-ray flash. 

Subject headings: galaxies: distances and redshifts - galaxies: high-redshift - galaxies: starburst — 
gamma rays: bursts 



1. INTRODUCTION 



The Swift satellite l|Gehrels et all 120041) promises 
the build-up of a significant sample of gamma-ray 
bursts (GRBs) with well-understood selection criteria 
usef ul for cosmological stu dies of high-rcdshift galax- 
ies JJakobsson et all l2005j) and the Hubble diagram 
l)Ghirlanda et al.ll2004afr The best way of securing GRB 
redshifts is from absorption line studies of the afterglow. 
However, for the first 56 Swift GRBs only 11 redshifts 
have been obtained. Alternatively, redshifts can be ob- 
tained from host galaxy spectroscopy although this re- 
quires that the host is sufficiently bright and has de- 
tectable emission lines. 

In this Letter we explore a third approach suitable for 
fainter host galaxies in case a spectroscopic absorption 
or emission line redshift has not been obtained. For a 
sample of 10 GRB host galaxies with spectroscopic red- 
shifts in the range 0.4 < z <2 and photometric measure- 
ments in more than 4 bands, [^hristcnscn et al. ( 2GG4a) 
found that photometric redshifts were consistent with the 
spectroscopic redshifts in all cases, within = 0.21. 
An advantage of the method is that it is independent of 
whether the host galaxy has emission lines. 

As the next step in validating the photometric red- 
shift approach we here predict a redshift from a host 
galaxy with multiband photometry which is sufficiently 
bright that a spectroscopic redshift can be measured, 
and thereby test the photometric redshift. Our target 
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is GRB 981226 which is currently one of the brightest 
hosts without a spectroscopic redshift. 

GRB 981226 was detected by the BeppoSAX satefiite 
on 1998 December 26.41 UT. It is consistent with being 
an X-ray fiash (XRF) with the fluences satisfying Sx > 
(Frontcract al. 2000). Despite intense optical and 
near-infrared (IR) follow-up observations i nitiated 6.5 
and 8.4 hours after the bu rst, respectively l)Kloselll99§: 
iCastro-Tirado et al] 119981) no optic al co unterpart was 
found fCalam a et 3.1.111991 lWozniaklll99a [Rhoads et alJ 
[1998; Bloom e talT 119981 ISchaefer et alJ 119981). The 
deep i?-band observations carried out bv Lind gren et alJ 
l(T999i) showed R > 23 mag at 9.9 hours after the GRB. 
This makes GRB 9 81226 a tvpical G RB without any de- 
tected afterglow (Tavlo r et alJ| ll99^ and close to b eing 
a dark burst jjakobsso n et alJ l2004: Rol et al.ll2005j) . 

Radio observations revealed a variable source 
at the position R.A.(J2000)=23''29"37!21, 
Dec(J2000)=-23°55'5 3'.'8 peaking ~ 10 days after 
the gamma-ray event ijFrail et alJ 119991) . Identification 
of the host galaxy was suggested based on the small 
angular separation between the radio afterglow and 
an extended object. High spatial resolution images 
from the HST/STIS show that the galaxy colours 
change notably over its s urface, with th e northern part 
being significantly bluer l)Hollandll2000() . The position 
of the radio after glow with resp ect to the galaxy is 
0'.'749±0'.'328 ( Blo om et all 120021) . which encompasses 
the blue northern part of the host. 

This Letter presents an analysis of all imaging data on 
the host available in public archives. Ground-based data 
in the BVRIJsKs bands and images from the HST makes 
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TABLE 1 
Log of the observations 



Instrument+filter 


date 


exp time (s) 




9nnn i o 




FORSI B 






FORSI V 


2001-05-16 


3x300 


FORSI R 


2001-06-19 


6x540 


FORSI / 


2001-05-16 


3x300 


FORSI / 


2001-08-13 


400 


FORSI / 


2001-08-19 


400 


FORSI / 


2001-09-22 


4x600 


ISAAC Js 


2001-09-22 


10x180 


ISAAC Ks 


2000-11-12 


30x120 


ISAAC Ks 


2001-06-07 


30x120 


STIS CL 


2000-07-03 


8265" 


STIS LP 


2000-07-06 


7909" 



''Total integration time. 

a multi-wavelength study appropriate to determine the 
host and hence the GRB photometric redshift. Some of 
the data have been presented previously in the literature 
in differe nt contexts. An ISAAC Ks band image was pre- 
sented inlLe Floc'h et a l. (2003), the HST/S TIS data in 
iHoUandl 1)2000(1 . and a deep R band image in lFrail et al.l 
the latter image is not analysed here. 

2. DATA ANALYSIS 

For consistency we re-analysed the photometry in all 
filters. The host galaxy was observed in the optical 
(BVRI) with VLT/FORSl and in the near-IR (Js and 
Ks) with ISAAC. The data were retrieved from ESO's 
public archive^. To take advantage of all available data 
on the host we also included images from the Hubble 
Space Telescope using STIS. The host was observed in 
July 2000, using a clear aperture (50CCD, or CL filter) 
and a long pass imaging filter, F28x50LP [LP filter)^. 
Tabled lists the dates of observations, number of inte- 
grations and exposure times. 

Data reduction was done using IRAF. All optical data 
were reduced using standard methods, i.e. bias subtrac- 
tion and dividing by a combined average flat-field frame 
obtained from twilight sky exposures. For the near-IR 
data reduction, sky subtraction was done by creating a 
sky image from the bracketing 8 frames from each indi- 
vidual night. After sky subtraction each image was di- 
vided by a normalised flat-field frame. Near-IR flat-field 
images were created by subtracting faintly illuminated 
flat-field frames from bright ones. Eight and ten of such 
images were combined and used for flat-fielding the Js 
and Ks images, respectively. Relative shifts were found 
using cross-correlation procedures before combining the 
individual exposures. 

Additional archive calibration files consisting of bias, 
sky flat-field, and standard star images from the differ- 
ent nights were retrieved. Photometric zero-points from 
ISAAC and FORSI were obtained by comparing the in- 
strumental magnitude of the standard star with tabu- 
lated values. All the zero-points calculated confirmed the 

^ http://archive.eso.org 

^ Images were obtained from 

http://www.ifa.au.dk/~hst/grb_hosts/intro.html (Holland 2001). 
ihe data are part ot the Cycle 9 programme GO-8640 "A Public 
Survey of Host Galaxies of Gamma-ray Bursts" . 



TABLE 2 
Host galaxy magnitudes 



Filter 


"iVoga(mag) 


niAB (mag) 


B 


25.74±0.37 


25.68±0.37 


V 


25.59±0.15 


25.63±0.15 


R 


24.81±0.06 


25.04±0.06 


I 


23.65±0.15 


24.09±0.15 


Js 


21.90±0.20 


22.84±0.20 


Ks 


21.10±0.27 


22.97±0.27 


CL 




24.79±0.10 


LP 




24.27±0.10 



Note. — Magnitudes of the host 
in the Vega and AB systems in vari- 
ous filters obtained with an aperture 
radius of 2". No correction for Galac- 
tic extinction is applied. 



values for the given dates that were found in the corre- 
sponding instrument web-pages. Extinction and colour 
terms are reported for each month on the FORSI web 
page and were assumed to be appropriate here. For the 
STIS data, zero-points were obtained from the STIS user 
manual. 

Aperture photometry was used to derive host photom- 
etry, and all magnitudes of the host listed in Table|21were 
derived using a 2" radial aperture. Corrections to larger 
apertures were found to be negligible. No optical after- 
glow was found for this burst and the observations were 
carried out 1.5-2.5 years after the burst. Therefore, no 
contamination of any significant level is expected for the 
photometry of the host. The magnitudes d erived here 
are consistent with R ~ 24.85 ± 0.06 mag in [ Frail et alJ 
(1999), Ks = 21.1 ± 0.2 mag inlLe Floc'h et all l|2003| r 
and CL = 25.04 ± 0.07 mag in'HoUandl l|2000() (using the 
same I'.'l radial aperture we find 25.00±0.05 mag). 

3. PHOTOMETRIC REDSHIFT 

We used the public p hotometric redshift code HyperZ 
llBolzonell a et al.*2 000(l along similar lines as in a series of 
papers on G RB h ost galaxies ( Gorosabel et al. 2003a ij; 
IChristeuseTi et abl f20fm Ic^rosahel et a,].l I200,^V As 
shown in these papers, photometric redshifts of GRB 
hosts can be determined to within Az — 0.21 when multi- 
band observations are available, and when spectral fea- 
tures such as the Balmer jump are bracketed by the ob- 
servations. 

Magnitude offsets between Vega and the AB system 
were calculated and added to the observed Vega magni- 
tudes listed in Table [3 The HST magnitudes are ob- 
tained directly in the AB system (see STIS Instrument 
Handbook), so the Table El fields corresponding to the 
HST Vega magnitudes are empty. To derive fluxes for 
the various pass bands the magnitudes were corrected for 
a Galactic extinction oiE{B — V) = . 022 m ag as derived 
from the dust maps of lSchlegel et"7nifT99l . Flux densi- 
ties in /iJy were calculated by f^, — io-0-4('"ab-23.9) j^j, 
an AB magnitude, toab, in each passband. 

Fluxes were compared to g alaxy template spec t ra cre - 
ated from the spectral atlas of lBruzual fc Charlotll|1993() . 
We used different galaxy templates o btained from a 
Salpeter initial mass function (IMF) ("SalDetei* '1953) 
(Sp55) and a Miller & Scalo IMF (M iller fc ScalQ.1979j) 
(MiSc79) with different star formation histories and a 
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Starburst at z.1 .1 1 , age=0.36 Gyr, Av.0.20 



TABLE 3 
Solar metallicity SED fits 




far 



HostOtGHB 981226 



1.0-10 1.5-10^ _ 2.0-10 2.5-10' 

observed wavelength (A) 

Fig. 1. — Best fit synthetic spectrum to the SED of the GRB 
981226 host galaxy using the ground based BV RIJsKs images. 
The template was created from a Sp55 IMF and a Calzetti ex- 
tinction law. The broad STIS filters are indicated by the squares, 
but are not included in the fit. Horizontal error bars indicate the 
widths of each band. 



characteristic timescale r, where r = corresponds to a 
starburst template, and r — > oo to an irregular galaxy 
template. Values of r in between these extremes cor- 
respond to various types of spiral galaxy and elliptical 
templates. These templates all assume solar metallicity. 
We also experimented using different extinction curves: a 
starburst extinction curve (Cal zetti et al. 2 000.), a Small 
Mage llanic Cloud (SMC) extinction curve l|Prevot et al.l 
I1984I). a Large Magellanic Cloud (LMC) extinction curve 
(iFitzgatrick 1986), and the Milky Way extinction curve 
(|Seatonlll97 9). 

In principle the STIS wide filters might be included to 
determine the photometric redshift. The CL filter ex- 
tends from 4000 to 9000 A and the LP from 5500 to 
9000 A, thus providing poor spectral information. Fur- 
thermore given the widths of the filters, the effective 
wavelength is sensitive to the assumed spectral template. 
Thus, as a first step, we decided to not include them in 
our analysis, and then discuss the impact that their in- 
clusion has on the Zphot determination. 

The best fit of the broad band SED was found with 
a starburst template at z = 1.11 with an age of 0.36 
Gyr and an intrinsic extinction of Ay = 0.20 mag as 
shown in Fig. ^ Using a Calzetti extinction curve and 
a Sp55 IMF gave the smallest reduced x^/d-o.f. as in- 
dicated in Table 131 Other combinations of template fits 
are also listed in Tabled In all cases a starburst tem- 
plate provided the best fit. From these fits we find Zphot 
= 1.11±0.06±0.10±0.21 (68%, 90%, and 99% confidence 
levels, respectively). In this case the Balmer jump is well 
sampled by the observations. 

The photometric redshift determination has a well de- 
fined minimum at z w 1.1 as shown in Fig. [3 The best fit 
values for the extinction, age, and redshifts are consistent 
independently of the extinction curve or t emplate used. 
Using the observed template spectra from iKinnev et al.l 
((19,9(J give consistent results for the photometric red- 
shift and extinction. Including also the HST bands in 
the SED fit, the resulting parameters (Ms, age, Ay) do 
not change. Specifically Zphot does not change at all, but 
the uncertainty for the photometric redshift decreases 
by ~50%. However, including the HST photometry in- 
creases the x^- A fit when including the LP band gives 
X^/d.o.f. = 1.315, while including the CL as well gives 
X^/d.o.f. 2.850. The reason for this is the calculation 



IMF ext. law Zphot {68%,99%) Ay (mag) x^/d-O-f- 



Sp55 

MiSc79 

Sp55 

MiSc79 

Sp55 

MiSc79 

Sp55 

MiSc79 



CalOO 

CalOO 

SMC 

SMC 

LMC 

LMC 

MW 

MW 



i.iilg: 
i.iil": 



06,-1-0.21 

04, -0.18 
06,-H0.2D 

06, -0.22 

07, -HO. 23 

03, -0.45 
06,-1-0.21 

05, -0.21 
07;-l-0.23 

04, -0.46 
06;-l-0.21 

06, -0.22 

07, -1-0.24 

04, -0.17 
06,-1-0.21 

05, -0.21 



0.20 

0.10 
0.12 
0.06 
0.16 
0.08 
0.14 
0.06 



0.399 
0.661 
0.427 
0.651 
0.426 
0.664 
0.432 
0.662 



Note. — Template SED fits using a Salpeter or Miller & 
Scalo IMFs, and different extinction curves. A redshift step of 
Az = 0.05 was used in our SED fits. Only ground-based photo- 
metric points were included in the fits. The derived properties 
in terms of best fit template type (Starburst), redshift, age (0.36 
Gyr), Mb = —20.25, and extinction are basically independent 
of the input parameters for the assumed IMF and extinction 
curve. A relative host luminosity L/L* fa 0.5 is derived assum- 
ing Mg = —21 mag. 




0.0 0.5 1.0 1.5 2.0 2.5 

redshift 

Fig. 2. — Reduced as a function of redshift with a minimum 
at z fs! 1.11 for the fit using the ground based data only. 

of the effective wavelength of the two HST/STIS bands. 
In Fig. |21the HST data points are shown by the squares. 

We investigated whether using templates with metal- 
licities different from solar values would have an im- 
pact on the output parameters. It is expected that 
the photometric redshift shows no significant difference 
while the well-known age-metallicity degeneracy would 
manifest itself. In Table 01 we list the results of fit- 
ting the observed SED with templates of metallicities 
of 0.2 and 0.4 solar, respectively. The template spectra 
were calculated using a Sp55 IMF with the GALAXEV 
code (fBruzual fc Charlotll2003j) . Generally the best fit 
templates have larger ages, while the photometric red- 
shifts are in good agreement with those in Table O 
Lower metallicities as well as templates constructed 
with a Chabricr (2003) IMF give consistent results for 
•Zphot- Since all the fits give similar values of x^j we 
can not disentangle this age-metallicity-extinction de- 
generacy using b road-band measurements alone (see also 
iBolzonella etal]l2000|) . 

To check for consiste ncy we used t he Bayesian pho- 
tometric redshift code fBe nited 12000(1 to estimate the 
photometric redshift. This code uses ernpirica l spec- 
tral templates for the fits l|Coleman et al.lTl98CI|) . Us- 
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TABLE 4 
Sub-solar metallicity SED fits. 



Z/Zq 


Zphot (68%,99%) 


Age (Gyr) 


Av (mag) 


X^/d.o.f. 


0.2 
0.4 


1 1 1 +0.09, +0.27 
-0.06,-0.22 
1 , 1 +0.06, +0.23 
-'-■-'--'--0.05,-0.21 


0.72 
0.51 


0.00 
0.16 


0.476 
0.399 



Note. — Here a Sp55 IMF and a Calzetti extinction law is 
used for creating the template. Only ground-based photometric 
points were used in the fits. 



ing the same ground-based photometric points as above 
we find 2phot,BPZ = 0.94 ± 0.13 (xVd.o.f. = 1.078), 
which is consistent with the results from HyperZ within 
1(7 uncertainties. Including the HST data points gives 
2phot,BPZ = 0.97 ±0.13 (xVd.o.f. = 2.457). The uncer- 
tainties reported here are Icr levels. 

4. DISCUSSION 

Based on multi-colour optical and near-IR photometry 
we have reported the first precise photometric redshift 
for a GRB host galaxy without a spectroscopic redshift. 
The host is sufficiently bright that the prospects of ob- 
taining a spectroscopic redshift, confirming or refuting 
our proposed value, appear feasible. We note however 
that we have been unable to detect emission lines in pub- 
lic VLT spectra in the ESO archive. In this connection it 
is interesting to note the relatively large mean age of the 
GRB 981226 host galaxy. This may indicate that emis- 
sion lines are not prominent in this galaxy and reminds 
us that GRB host samples with spectroscopic redshifts 
may be biased towards emission line galaxies. The VLT 
spectra were obtained at a relatively low resolution (>10 
A), and higher resolution data should be obtained to test 
the photometric redshift. At z « 1 the optical emission 
lines could well be contaminated by residuals from sub- 
traction of the skylines. Instruments such as ESI, MIKE, 
or the next generation instrument X-shooter could reveal 
such cases. 

In addition to the photometric redshift, the SED 
also allows us to infer physical properties of the host 
galaxy. Assuming a redshift z = 1.11 and the cur- 
rently favored flat cosmology with = 0.3 and Hq — 
70 km s~^ Mpc^^, the luminosity distance is 2.3 x 
10^^ cm. By interpolating a power law function between 
the observed B and R band AB magnitudes (corrected 
for extinction), we find that the rest frame 2800 A Ultra 
Violet (UV) flux is 0.27±0.06 ^Jy. This corresponds to a 
star-formation rate (SFR) of 1.2±0.3 Mq yr~^ using the 

3 BATSE detected a short GRB (T9o=1.7 s, BATSE trigger num- 
ber 7281) which has been associated with GRB 981226, but the 
BATSE burst was triggered one hour after the BoppoSAX detec- 
tion. Hence, these two bursts are different, and the peak energy 
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conversion between UV flux and a global SFR (jKennicutti 
11998ft . The absolute magnitude of the host galaxy is 
Mb = —20.25 mag, which implies a specific SFR of 
2.4±0.6 Mq yr-i (L/L*)-i, where L* corresponds to 
the luminosity of a galaxy with Mg = —21 mag. This 
specific star formation rate is smaller that the UV based 
SFRs in a sample of 10 GRB hosts ijChristensen et all 
I2004a|) . As for other GRBs there is no indication that 
the host galaxy is strongly affected by extinction. 

GRB 981226 had the interesting properties of 
being consistent with being both an almost-dark 
GRB llJakobsson et all 120041) and an X-ray flash 
l)Fronter;^talTl2'00nl) '^. We can use the inferred Zphot 
to estimate the peak energy for the BeppoSAX burst. 
Integrating a Band function and adding four power law 
functions to represent the X-ray to 7-ray spectrum at 
different time segme nts with the parameters given in 
iFrontera et al.l l|2000() . the total time integrated fluence 
is 1.9 X 10~^ erg cm~^ and at z = 1.1 1 the isotropic 
energy release is Eiso = 5.9 x 10^^ erg. IFrontera et alJ 
l)2000fl divided the light curve into 5 segments, one of 
which had £^pcak— 61±15 keV, the others less than 
10 keV. The Amati relation (iGhirlanda et alJ200'4aj) pre- 
dicts £;pcak = 77.6 keV and E°^ly^ = 36.9 keV, consis- 
tent within 2a with the BeppoSAX observations of the 
peak energy. With small uncertainties in the luminos- 
ity distances, accurate photometric redshifts for other 
GRB hosts will help to constrain cosmological param- 
eters through the Ghirlanda relation. The photometric 
redshift is higher than those reported for other XRFs but 
still lower than the median redshifts for GRBs, consistent 
with the hypothesis that XRFs may be off-axis GRBs or 
dirty fireballs which both predict lower mean redshifts. 
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